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The problem of surface effects at a fluid/force field boundary is investigated. A classical simple fluid with a locally 

^^ ' introduced field simulating a solid is considered. For the case of a hard-core field, rigid, exponential, realistic, and 

macroscopically smooth boundaries are examined. 

O ' Two approaches to this problem are analyzed, which have been developed independently from each other over a 

U ' long period of time. With some degree of arbitrariness, they can be referred to as "adsorption" vs "surface tension" 

II , or "cluster expansion" vs "pressure tensor" . 

The "adsorption" approach is used to obtain a series in powers of the activity for the specific surface omega 
potential (grand potential) - 7. This expression is similar to the Mayer expansion, but, unlike in the case of pressure, 
the integrals of the Ursell functions contain factors which depend on the particle/wall interaction potential. These 
tS^ ' factors are written in the general form for an arbitrary term of the series, and have a fairly simple structure. In the 

yf , case of a hard wall, the coefficients of the series reduce to the first moments of the Ursell functions taken over certain 

regions. Equivalence of the "adsorption" and the "local pressure" approaches is established. 

The "surface tension" approach is used to expand the Kirkwood-Buff formula for the surface omega potential of 
the field/fluid system to the arbitrary localization of the dividing surface. The surface terms identification depends 
C^ ' on the particular form of the fluid particle/solid interaction potential, and on the formulation of the problem. As 

ty^ ' a function of an activity, the surface terms can begin both with linear term and quadratic one, depending on the 

situation, and this corresponds to the presence or absence of adsorption in the ordinary sense. "The surface tension 
C^ ' coefficient" breaks up into the term proportional to the Henry constant, depending on the dividing surface position, 

and universal nonlinear surface coefficient. 

Using methods of statistical mechanics, it is shown that the derivative of the tangential component of the pressure 
' "sJ ' tensor with respect to the chemical potential coincides with the near-surface number density on average over the 

H , transition region, that has two consequences. 

Firstly, it proves complete identity between "tension" and "adsorption" approaches in the domain of their existence. 
Thus, the results of one line of investigation can be used to develop the other. In particular, the "adsorption" method 
obtains the expression for the surface terms in quadratures, and the "tension" technique - the analog of the Mayer 
expansion. 

^ I Secondly, it gives the near-surface virial expansion, which determines the exact equation of state of near boundary 

^\^ . "two-dimensional" fluid. The tangential component of the pressure tensor averaged over the transition region plays 

jir ' the role of pressure, and the average number density - the role of number density. 

^ -, . . . 

Q 1. Introduction also corrected fS','^'. 

^^ , A significant drawback of this approach [SUB] is that the 

^~~^ ' The surface terms at the solid/fluid interface for the case computation algorithm is extremely complex or, in other 

J> ,; of "gas adsorption on a solid surface" have been studied words, there is no universal expression for the general 

in terms of the statistical mechanics for a fairly long time term of the series. 

[mi]. Most papers consider a solid as a constant field, as The second problem is related to the Henry constant, 

C^ ' the first step of a more complex problem. Here we will and here we observe three cases. 

not deal with the problems related to the motion of solid In the most exotic one, the interface locates at the first 

particles. layer of solid atoms [5J|ll|Slin] . 

More correct calculation [S], as shown later, corre- In some studies this constant is strictly zero, while 

sponds to a hard solid. Attempts were made to combine in many other papers, it is proportional to the excluded 

these results [5]. volume pi lTUlfTT] . 

In subsequent studies, the initial consideration [5] was The equations used in some papers 0, [12] give nega- 

extended to the case of an arbitrary particle/surface inter- tive values for the number of surface particles under cer- 

action potential [7HS]. Inaccuracies of early studies [311] tain conditions. 

related to the integration domains of the series terms were These problems are due to the fact that the technique 
of arbitrary localization of the dividing surface has not 
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This question is closely related to the issue of interpre- 
tation of the surface number density, so, it is also not fully 
resolved. 

Another reference point in the interface problem is the 
Kirkwood-BufF formula [13] for the surface tension of the 
liquid/gas system 
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[P - P\x)] dx. 



(1) 



Here P is the pressure in the homogeneous phases, P' 
is interpreted as the tangential component of the pres- 
sure tensor |14[|15) , and x is the coordinate normal to the 
interface. 

In addition to the force-based consideration [T5] . 
a similar formula was derived earlier using a quasi- 
thermodynamic analysis [TB]. In this case P' has the 
meaning of the equilibrium pressure corresponding to the 
number density in the transition region. 

Equations such as ([T]) are also applied to the solid 
(field)/fiuid interface [TTllTTllTS] . (In some cases, the nor- 
mal component of the pressure tensor - P„ was extracted 
from P |11[[T71[T^ . In terms of further consideration, this 
manipulation seems meaningless.) 

An application of this equation to the field/fluid inter- 
face, in our opinion, requires further validation. 

Particularly the question of integration limits in equa- 
tion ([1]) is not quite clear. In the case of liquid/gas sys- 
tems the integration is performed from — oo to +cxd, how- 
ever, it is impossible for solid (field)/fluid systems due 
to P' zeroing inside the solid phase. For some reasons, 
possible variants of generalization [20] have not received 
logical continuation. Thus, here we also have the above- 
mentioned problem of the dividing surface localization. 

The "mechanical" consideration (section [3. lOp has one 
significant drawback: it deals with the tangential pressure 
which is distorted by the transverse wall (section 14. 3p . 
This issue is not adequately covered in the literature [17] . 

Finally, from the viewpoint of analytical space expan- 
sion, and the improvement of the problem "ecology" , it 
is of interest to establish whether the "adsorption" di- 
rection [THSIIHIIS] is equivalent to "surface tension" direc- 
tion [I3lll5lll7lll8ll21) concerning force field/fluid interface. 

Some results, related to this contribution were obtained 
in our previous work f22L In particular, this regards to the 
technique of partial localization factors, which generalize 
the notions of Boltzmann and Urcell factors. 

2. Primary definitions 

2.1. Canonical ensemble 

The probability density of finding a given spatial config- 
uration of a particular set of particles [U p. 181] is given 
by the expression 



^l...k — 
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exp{-l3Ui ^)drk+ 



.dr 



N- 



(2) 



Here N is the number of particles in the system, /3 = 
l/ksT^ ks is the Boltzmann constant, T is the tempera- 
ture, U^ pf is the energy of interaction between the parti- 
cles, and V is the volume of the system. The integration 



is performed over the coordinates of the particles of the 
ensemble Vk+i, . . . , r^. Zn is the configuration integral 



Zn = exp{-(3U^ j^)dri...drN. 



(3) 



Going over to the distribution functions for an arbitrary 
set of particles, we obtain 
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where g^j^ j, is the probability density of finding a given 
configuration of k arbitrary particles for the canonical en- 
semble. 

2.2. Grand canonical ensemble (GCE) 

Let us average equation ^ over fluctuations in the num- 
ber of particles by applying to both of its sides the oper- 
ation X]?/=o-^JV' where 



z Zn 



pV 



(5) 



is the probability for the GCE to have a definite number 
of particles N inside the volume V. Here z is the activity: 



~A^' 



(6) 



where /i is the chemical potential, A = /i/v^27rmfcsT, h is 
Planck's constant, m is the particle mass, and Sv is the 
large partition function of a system of volume V: 
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We have 
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X / exp(-/3L/]^^]^_^j,)drfc+i...drfe+Ar|. 



(Tq \ J, specify the probability density of finding a cer- 
tain configuration of k arbitrary particles in the GCE. For 
an ideal gas, (Iq^ i^ — q'', where g = N/V is the number 
density. 

2.3. Ursell factors and partial localization factors 

(k) 

The Ursell factors U^ \, , which are also called cluster func- 
tions, appear in the well-known expansion of the pressure 
in powers of the activity [4, p. 129], [221 p.232] 



oo ^ ^ 

P(z, T) = zkeT + keTY,^ / ul^']^dr2...drk. 

fc=2 ' '' 



(10) 



They are also contained in the corresponding expansion 
of the number density 
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which is obvious from the relation |j 

dP 
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The Ursell factors decay rapidly as any group of parti- 
cles, including a single particle, moves away. 

Because ([TO)) is, in fact, a Taylor series expansion and 
in view of the local nature of the Ursell factors, we can 
write the equation 
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which is valid on macroscopic scales, and where the 
derivatives are taken at the point z = 0. Here 6{r) is 
the Dirac delta function. Naturally, in this case, we as- 
sume that the index k does not reach macroscopic values. 

Some other properties of the Ursell factors are described 
in Appendix lA.il 

The partial localization factors Sj™j^_^j, are the hybrids 
of the Boltzmann and Ursell factors. They are contained 
in the expansion of the distribution functions in powers 
of the activity 
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are characteristic functions. 

(m^k) 



The properties of B\^ 



.m+k 



are described in Appendix 



2.4. Presence of an external field 

The configuration integral of the inhomogeneous closed 
system is given by the expression 

N 



which is obviously the large partition function of the sys- 
tem in the presence of an external field. 

Using the properties of the fractional generating func- 
tion dni) (see Appendices IA.2I [13. 1|) . we obtain the analog 
of (HID 

cxD „ „ n+1 

e(ri,z) = Z01 + Z01 ^ ^ / [ n ^»]z^i':!.+Wr2...dr„+i, 

(18) 
where £'(ri, z) is the number density in the presence of an 
external field and 



9i = exp(-/3wi) = exp[-^u(ri)] 
and the analog of (IT4l) 
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3. System in a force field 

We consider a statistical system in a field whose range is 
much smaller than the size of the system and which is far 
from the boundaries of the system. That is, the case in 
point is rather "the field in the system" . 

We also assume that the potential of the fields consid- 
ered increases/decays rapidly enough to ensure the con- 
vergence of the corresponding integrals for the motion in- 
ward and outward the solid, and the interparticle poten- 
tial decays rapidly enough to provide the convergence of 
zero and first moments of the Ursell factors (i.e., we do 
not consider the Coulomb potentials in this paper). 

In some cases, we assume that the range of the solid 
field has a finite radius, so we speak about a transition 
region between the solid and the fiuid. This approach is 
of no fundamental nature and will be used to simplify the 
consideration. 

Finally, unless otherwise stated, we assume that the 
region of the transition from the solid to the fluid is much 
smaller than the linear sizes of the region of the field, 
and that the field profile along its gradient is the same 
everywhere. 



Zn = j exp(-/3^ w, - /3U^^)dri...drN, (16) 3.1. Basic equations 



where Vi is the energy of the interaction of the i-th particle 
with the field. 

For the GCE, we introduce the quantity 



Taking the logarithm of the ratio of the partition functions 
([T71) and © gives 
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^To simplify the formulae, we do not indicate temperature con- 
stancy; it will be assumed in all derivatives. In addition, differ- 
entiation parameters will not be indicated in cases where this is 
obvious, for example, when they are set by the opposite side of 
the equation. 



expi-P^Vi) - 1 Ul%dri...drk 



where the logarithmic form of the generating function for 
the Ursell factors is used. 



Under the above assumptions, we can obviously omit 
the symbol V in the integrals over the coordinates of par- 
ticles in (I^TI) and consider the integration to be performed 
over the infinite space. 

We rearrange (PTjl as follows: 



lnS^-lnSv = -^ 
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ip{ri) = 1 - exp{-(3v,) = 1 



(23) 



As can be seen from (fTU|) . the first sum in (P^ is pro- 
portional to the Mayer series for the pressure of a homo- 
geneous system. 

The second term on the right side of (0^ is proportional 
to the area of the surface which separates the field region 
from the rest of the system. Indeed, its structure is such 
that at least two particles are always on the opposite sides 
of the "field boundary" . Expanding the product in any 
term of it, we have 
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where "2 < j < t. With the first particle fixed, the inte- 
grals of the type of ((M)) are determined by a local region 
near this particle because of the locality of the Ursell fac- 
tors. Consequently, J ...dr2---drt do not depend on the 
displacement of the first particle along the boundary of 
the system. 

As the first particle moves from the boundary into the 
field region, the integrand of ([M|) decays rapidly due to 
the Boltzmann factor 9i . 

As the first particle moves from the boundary into the 
fluid region, it decays rapidly because of the fixing factors 
(fi and the locality of the Ursell factors. We can say that 
the Ursell factor "glues" functions 9i and ipi, i.e., the fluid 
and field regions. 

Integrating over the surface, factoring the area out of 
the sum, and after simple manipulations, we obtain 
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-PV + P ifdr + vA, 
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where Vt^ is the omega potential (grand potential) of the 
system with the field embedded in it, A is the area which 
bounds the field region, and we introduce the notation 
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(26) 
where x\ is the coordinate directed along the field gradient 
and it is chosen so as to satisfy the condition dx\ > 0. 

Recall that the integration is performed over the infinite 
space. 



The first terms of series, which is analogous to ([26|. 
were calculated in the set of papers [HIIS] based on a topo- 
logical approach [5] extended beyond the potential of a 
hard wall. 

Since the Ursell factors depend only on the relative con- 
figuration of particles, in (j26p we can perform integration 
over dxi. Making the change of variables ri — ri,r^ — 
ri — ri,i — 2, . . . k and using the invariance oiUl \. under 
translations, we get 



iy{z,T) = ksT 
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(28) 



is the function of the variables X2, ■■■Xk, which is symmet- 
ric under permutations of particles, and has the dimension 
of length. Obviously, in contrast to Z//-[ j^, it is not invari- 
ant with respect to translations since it depends on the 
location of the external field. 

Note that in ((27)1 the variables are separated: the ex- 
ternal potential is contained only in /C^^^) and the inter- 
parlicle one - only in lyl^''\ 

Expressions ([25)1 and (|?7)l are basic for the further con- 
sideration. The first term on the right side of ([^5)1 has 
purely bulk properties, and the third one has purely sur- 
face properties. The second term is of a mixed nature and 
has both bulk and surface properties, as is easily seen from 
its structure. 

It is obvious from the foregoing that the attempts to 
represent the f2-potential as a sum of uniquely defined 
bulk and surface terms are, generally speaking, incorrect. 
There should be a free parameter - the position of the 
dividing surface. 

The function ly plays an important role in the problem 
considered. Other representations of this quantity will be 
given below. Formula ([27| is, in fact, an analog of the 
Mayer expansion ()10p . which, generally speaking, cannot 
be further simplified. Note that, as is easily seen, the 
expansion of the function v{z^T) in a series of z begins 
with a quadratic term. Thus, it is logical to call this value 
a nonlinear surface coefficient. 

3.2. Variants of accounting for surface effects 

In the derivation of (|25)l - ()28ll . we assumed that the po- 
tential of the external field is independent of the displace- 
ment along the boundary. Suppose, moreover, that it has 
a hard core, and, for simplicity, in some cases we assume 
that the transition region is strongly restricted to a cer- 
tain length. 

We rearrange ([25)1 as 

x' CO 

Vt = -P{V - V) -PA I 0dx + PA I (pdx + vA, (29) 

— CXj x' 

where x' is some arbitrary point within or near the tran- 
sition layer that determines the volume V' , and x is the 



coordinate directed along the field gradient. (With the 
same result, the integration region in equation (j29p can 
be restricted to the transition layer because outside it the 
integrands vanish.) 

We see that there is some degree of arbitrariness in the 
differentiation between bulk and surface terms. Actually, 
this arbitrariness exists only in the transition layer or near 
it. Indeed, it is clear that otherwise we have the banal 
compensation of two identical terms of opposite signs. 

We denote the volume of the hard core by Vi , and the 
volume of the transition region by Vt- As noted above, 
Vt<^Vi. 

From (P5)) we obtain at least three variants of account- 
ing for surface effects for different values of x' : 



n 



-P{V -Vi- Vt) -PA I 



9dx + uA, 



(30) 



3.3. Surface number density 

Differentiating (|29p with respect to the chemical potential, 
in view of 

and using ([T^. we obtain 



N = q[V~V') + qA / Odx-gA / (pdx-A(3z-^, (37) 



N^Nb + N„ (38) 

where iVf, and TV, are the number of bulk and surface 
particles, respectively. For them, we have 



m = g{v - r), 



(39) 



n = -P{V - Vi) + PA / ifdx + vA 



I' 



(31) 



N.^gA I Odx - qA I ipdx - A/Sz-^. (40) 



n 



-P{V-VQ) + iyA, 



(32) 



where It is the length of the transition region and Vq is 
the volume of the solid bounded by the surface defined by 
the condition 



Odx = / (pdx 



(33) 



On the left side of this equality, the integration is per- 
formed outward from the field region, and on the right 
side, inward to the fluid, and it determines the position 
of the surface in question xq- 

Variant (j30p corresponds to x' lying on the boundary 
of a homogeneous fluid, i.e., it assigns surface effects to 
the solid, while variant ([51]) corresponds to x' lying on 
the boundary of a homogeneous solid and assigns surface 
effects to the fluid. Variant p2p sets the linear surface 
terms equal to zero, using the property ((33)) . In this case, 
x' can be both within and outside the transition region, 
depending on the form of the external potential. 

In computational studies, one often uses the fourth vari- 
ant in which the integration is extended to the inacces- 
sible volume near the surface to keep the volume of the 
entire system. This example will be discussed below. Of 
course, the other variants based on formula ([2^ can also 
be considered. 

The results obtained above will be discussed below, and 
here we present the traditional expression for (P^)) 



n = -P{V-V')+-f{x')A, 



(34) 



7(x') = -P / edx + P ipdx + iy (35) 



where 



is the general form of the specific surface fi- potential. In 
this expression, the term linear in pressure is separated 
from the nonlinear one - v. 



Going to the surface number density 



Qs 



A 



and taking into account that ([T5|) leads to 



diy 



Pz-g^ = I [g{x) ~eg]dx, 



after simple manipulations, we finally obtain 



(41) 



(42) 



g{x)dx + / [g{x) — g] dx, 



(43) 



where g{x) is the number density near the field boundary 
as a function of the coordinate directed along the field 
gradient. 

As will be seen from the foregoing (|43p agrees with or- 
dinary equation 

--!■ <«' 

Formula (|15)) generalizes the expression for the surface 
number density [TJIHIIII] to the case of arbitrary location 
of the dividing surface. 

Thus, the contribution of the solid phase to the surface 
number density is defined by g{x), and that of the fluid - 
by g{x) - g. 

It is easy to see that actual contribution to (H5)) comes 
from two regions: the region of the interaction between 
particles and the surface and the adjacent region of near 
surface distortions (oscillations) of the number density. 
(This two regions can intersect.) 

For the two cases for a hard solid considered above (sec- 
tion |33)j we have 



= / [q{x) - g] dx 



(45) 



and 



oo 

Qs = -Rg+ / [g{x) - g] dx, 



(46) 



Now return to basic relations. It follows from (|29| that 
n^-P{V -V')-PAKHix') + :^A, (55) 



where R is the particle radius and the point xq is placed 
at the coordinate origin. 

Expression (|43p has a simple meaning. Relying on the 
obvious relation 

(47) 



(48) 
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-PKh{x) + V. 



(56) 



N = I g{r)dr 

V 

and using (|55|) and (|39[). we obtain 

TVs = AT - iVb = / ei(r)dr - g[V - V"') 



This equality is one of the most significant expressions 
in the case under study. It combines the three most 
important surface parameters - the specific surface $7- 
potential 7, the Henry adsorption constant Kh, the non- 
linear surface coefficient v, and one bulk parameter - pres- 
sure. 

Differentiating ([55)) with respect to the chemical poten- 
tial and using (|35]), we obtain 



-|-oo 



g{r)dr + / [g{r) - g] dr, 
v v-v 



gs = gKnix') + / [g{x) - 9g] dx, 



which, under the assumptions made, is identical to (j43p . 
Thus, this expression is valid for all fluid densities. Cor- 
respondence of (P^l) and ([35]) demonstrates the internal 
consistency of this approach. 

3.4. The Henry adsorption constant 

For the Henry constant 

,. Qs 



Km — lim 



e^O g 



(49) 



(57) 



which is another expression for the surface number density 
for any position of the dividing surface. Here, the parts 
linear and nonlinear in the bulk density are clearly sepa- 
rated. This expression is more convenient than its analog, 
derived on the topological basis for a particular value of 
x' and a particular kind of particle- wall potentials [9] . 

Thus, we can consider a partition of the surface quan- 
tities on the basis of linearity with the same rights, as on 
the basis of the localization. 

Indeed, if we define a linear surface number density 



([101) gives 



K 



ex.p{—/3v)dx — / [1 — ex.p{—/3v)]dx. (50) 



H 



This expression generalizes the earlier obtained one [5] 
to the case of arbitrary localization of the dividing surface 
and an arbitrary law of the increasing potential when we 
move inward the body. 

The Kh dependence on x' is universal and can always 
be extracted explicitly. Indeed, subtracting one value of 
this constant from another we ascertain that 



gsj = gKnix') 
and nonlinear one 

-I-CX3 

Qs,n = / [gix) - 9g] dx, 



Kh{x')^x' -x" + Kh{x"). 



(58) 



(59) 



then, taking into account ([T2|) and \A2\ . we get ordinary 
relations 

„ - - ^[^^-(-')] (60) 

(61) 



(51) 



Choosing the value x" — xq defined by ([55)1 , and taking 
into account that 

Kh{xo) = 0, (52) 

we obtain 

Kh{x')^x' -xa{T), (53) 

and the problem of Kh calculation reduces to the deter- 
mination of Xq . 

From (|33p . it follows that 



IJs,l — 


9/i 


and 


du 


Qs,n = 


dfi' 


in agreement with ([44)). ([56]). 




Naturally 




Qs = gs,i 


+ Qs.n, 



Xo 



{x)dx + / tp{x)dx 



(54) 



(62) 

and ([55)) - ([5^ define the surface number density in the 
most general form. 

3.5. Abruptly changing field 

Let us return to the general equation ()25p and consider 
a wall with a hard solid potential. (But we can still con- 
sider the case of an arbitrary interparticle potential.) The 
potential energy of particles in the field of the wall is ex- 
pressed as 



Surface passing through the point xq by virtue of ([52 
can be called zero adsorption surface. 



w(n) 



+CX) (n e Vq) 
(n ^ Vo) 



(63) 



In this case, the functions ipi ([23]) take the form of the Using the change of variables rj = r2,r[ = r^ — r2, for 
characteristic functions (J15p i = 3, ...t and omitting the primes, we arrive at 



t p 

{z,T) = kBTj2 777^ / X2ll^o%..tdr2-.drt, 

*_o I^V^ ~ ^J- Jxo.>0....x,>0 



^^ ^ ^r- (64) „r. T) = k^T\- " ../ii*) 



(^2 ^(^ ~ ^-'^ "'2;2>0,...a;t>0 

Thus, the sohd occupies a macroscopic volume Vq inside (72) 

V. As noted above, we assume that the solid is far from since the Jacobian of this transformation is unity, 
the boundaries of the volume V. In this case we used the invariance of W^*) under parti- 

Equation ([25)) leads to cles permutations and under spatial inversion. 

For the Henry constant ([50]) in accordance with ([55]) in 
^" ^ -{V -VQ)P{z,T) + Aiy{z,T) = nv-Vo, (65) this case we obtain 

where for v{z,T) before the integration over xi, we have Kh — x — Xstep, (73) 

the expression 

where Xstep is the step location, which generalizes case 



v{z, T) = kBT2_^ —r (66) In the above, we considered the case where the volume 

t=2 ' of the solid is limited by condition ([63]) . i.e., in fact, we 

/P * 1 r -1 increased the volume of the field by the quantity AR. 

xY" [l - n(^ ~ ^^^Y"^Ul*'^dxldr2...drt, (This corresponds to condition §^ for the field of a hard 

*=2 solid.) In some cases, it is reasonable not to include this 

, . ^1 1 ■ i T 1 i xi r region in the solid. Then, we obtain the expression 
where xi is the coordinate perpendicular to the surface 

given by the boundary V^r^ and n^ ^-{V -Vq + AR)P{z, T) + A[RP{z, T) + v(z, T)], 

. (74) 

Vo _ , /^o _ J ^ ('"i G ^o) , s in which the bulk part of the solid coincides with the 

* * I 1 (ri ^ Vb). field region, and the surface term has a part linear in 

number density |ll[Tni[II]. Naturally, this corresponds to 

The function f/'i^ localizes the i-th particle in a volume the Henry constant ([73|) x' - Xstep = ~R- 

V, and the function x^ localizes this particle outside this In conclusion, a few words about hard spheres, 

volume. A system of hard spheres at a hard wall is an object 

The first terms of the expression which is analogous to of a constant attention as a convenient model exercise. 

([M)) were obtained using a diagram technique [Sj. The For the second surface virial coefficient calculations give 

series as a whole was examined in the construction of the minus ttD'^/S, where D is the sphere diameter, and for 

open statistical ensemble [2211111 ■ *he third one - minus Ugn'^Dy 1680 ^\7M\- 

We perform integration over xi in ly. The functions 9 The calculations, carried out in the present work on the 

in ([281) also take the form of characteristic functions basis of ([Ml) , dZD) and independently ([HH) show complete 

coincidence of results with these data. 

fl ^ _ , , / 1 (:r, < 0) 

'^t ^- Xt ~ X\^i) — <, „ (x > Q) ^'^- Macroscopically smooth field 

We assume here that the characteristic size at which the 
and the product 6i in the integrand in ^ is obviously field undergoes significant changes is macroscopic. Sub- 
determined in this case by the maximum Xi stitution of ^ in ^ and integration over r2...rfc give 



+ 00 



/i'rfe'^^ f Xixi)[l-x{xi+ma^{x2...xk})]dxi. (69) K^'^) = E fi g^ / {^(^) - [^(^)]'} '^^ (^5) 



^k Qkp 

k=2 



Thus, we obtain 

/i';rfe'^=max{0,a;2...Xfc}, (70) 

and for v{z,T) we have the expression 



(z,r)= / [ep{z)-p{ez)]dx, (76) 



.(z,r) = fcBr^^^^-^ (71) resuh 



where the terms linear in the activity are added to both 
integrands. 

Substitution of this expression in ([25]) yields the logical 



X /' .2<..A....^r., n- = - J PiOz)dr = - J Pize-^^dr = 

Jx2>{0,X3...Xt} V V 

where the symmetry of the Ursell factors under particle — / P[^ — u{r)\dr — — P{r)dr, (77) 

permutations is used. ^ ^ 



where we used the well-known expression for the chemical 
potential in the field. 

From (|50l) . ([5E| we obtain 



X oo 



(78) 



Expressions ([75|. (|75)) may be useful in studies 
of surfaces with macroscopic thickness or in quasi- 
thermodynamics. The typical structure appearing in (I78p . 
as will be seen from the following, is universal for 7. 

3.7. Exponential field 

Consider the "wall" -fluid potential of the form 



v{x) = t;oexp(--), 



(79) 



which used in the numerical calculations |10) , and deduce 
the surface functions for it. 

Using (|50p and the exponential integral expansion 



Ei(£) = C-t-ln(~e) + ^— , [e<0] 






(80) 



some distance far enough from the interface in the region 
x < and for the rest is identical to ([7^ . This distance 
is determined by 9 proximity to zero. Such a substitution 
does not affect the calculation, but will comply with the 
requirements. 

Taking into account the relative simplicity of the ob- 
tained relations it is not improbable that the "wall" with 
the present potential may be used along with the hard 
wall as a convenient model system. 

3.8. Microinhomogeneous field 

In this section we consider briefly the situation when the 
field depends on the displacement along the surface. We 
assume that the simulated solid is crystalline, and thus 
the field has a periodic surface structure. 

Considering the integrals over the coordinates of the 
first particle that are parallel to the surface, we arrive at 
their averaging over the surface unit cell. 

Thus, Lp in ((25)) is replaced by 



M = 



1 



O-cell J cell 

and /i \ in (P7|) is replaced by 



(fidydz, 



(86) 



where C ~ 0.577 is the Euler constant, in accordance with 
([55)1 . for the Henry constant we obtain 



+00 



fik-l)\ 



un") = / dx. 



Kh = x' - l{\n 



knT 



-C). 



(fe-1) 



In a similar way, for the function f^ j. 



f(fe-i) 

J2...k 



/In 



1 + Y^eM-^) 



(81) 
we get 

(82) 



1 



Icell J cell 

+00 



rAl-\[e{ri+ri)\dyidzi (87) 



1=2 
k 



I {e{r,)) - {e{r,)\{e{r,+n))dxu 



When the rn-th particle moves far away from the inter- 
face inward to the field region, Xm < and 



Ak-i) 

J2...k 



(83) 



which does not disturb the convergence of the integrals in 
(|27)) due to a stronger decay of the Urcell factors U[ \.. 

When / -^ expression (15^ reduces to ()70p . that is 
easily seen from the behavior of the exponents in (|5^ . 

If parameter I takes a macroscopic value, then, using 
(|T^ . we obtain 



f2...k^^o%..kd'r2-drk = ipink—^ 



and taking into account ()27p . we arrive at 

z'^Xnk d^P 



v{z,T) = lY,- 



(84) 



(85) 



1=2 



where aceii is the area of the surface cell, and an averaging 
over y,z at a given depth of x is indicated by triangular 
brackets. 

Equations (|29)) - (|33)) retain their form, with the re- 
placements of (/? — ^ {(fi) and /2 ^ ~^ (/2 ~k )■ Expres- 
sions (1531) . ()56p seem to retain their original appearance. 

In the case where the contact surface has crystal faces 
which are differently oriented with respect to the crystal 
axes, the surface terms must be summed over these faces. 

3.9. "Local pressure" 

Integration of (|55)) with respect to the chemical potential 
yields 



rj = - 



Nd^Ji', 



k=2 



naturally, in accordance with ()75p too. 

Despite potential ()79p is beyond the initial setting of the 
problem (the hard core presence), the above results are 
valid at least in some temperature region. To prove this 
conclusion, we consider a potential that becomes -|-cxd at 



where we leave out an arbitrary function of volume, tem- 
perature, and field which arises from the integration, since 
in the limit of low densities the J7-potential must tend to 
zero. Using (jTT)) . we arrive at the relation 



n^- f P*{r)dr, 



(89) 



where 

P*{r) = I g{r)dti', 

— oo 

which, in this case, is equivalent to 



g{r) 



dP*{r) 
9/i 



(90) 



(91) 



We see that the number density for the system in the 
field and P* are linked by exactly the same relation as 
the usual number density and pressure for a homogeneous 
system p^. 

A quantity equivalent to the P* at a small number den- 
sity was introduced when fluctuations in inhomogeneous 
medium were calculated 20 and we shall also call it a 
"local pressure" . 

Under the conditions of the Mayer- type expansion va- 
lidity, relation ([T5| implies 

P*{ri,z)^kBTzei (92) 

which also resembles the standard expansion of pressure 
in powers of the activity ([TUl) . 

It should be noted, however, that P*{r) and g(r) do 
not correspond to one another in a macroscopic sense, 
i.e., they are not linked by a quasi-thermodynamic re- 
lation. It is clear that this would be valid for macro- 
scopically smooth fields, for which one can apply ([T^ . 
In this case, integration over the coordinates yields the 
pair P[n — u{r)] and q[^ — u{r)], related by an ordinary 
equation of state. Otherwise, the functions 9 distort the 
integrals of the Ursell factors in (|18p and ([M|) . and we 
have a different functional relationship. 

This consideration, in particular, shows the range of 
application of the quasi-thcrmodynamic approach [16] . 

It is seen from ((90)) that P* rapidly tends to P when 
the point of observation moves away from the field region. 

We emphasize that equation (|89l) is absolutely general 
and valid for a system in a force field of arbitrary config- 
uration. In the case of macroscopically smooth fields, it 
transforms into equation (1771) . 

Similarly to the above consideration of the case with 
the number density, it is possible to reduce equation ([M|) 
to 

n^-P{V-V')- f P*dr+ I {P-P*)dr (93) 
v v-v 

by simple manipulations. 

Integration over the surface yields 



n 



X oo 

-P{V -V')-A I P*{x)dx + Af[P~ P*{x)]dx. 



(94) 
Just as it should be, the derivative of (IM|) with re- 
spect to the chemical potential gives the expressions for 
the number of particles: bulk ([5^ and surface (P5)) ones. 



For the specific surface J7-potential 7 we obtain 



7(x') = - / P*{x)dx + [P- P*ix)]dx. (95) 



An equation analogous to (j95p was also obtained when 
considering fluctuations in piped system pO] . 

It is obvious that ([M]) ensures that p4p is valid. 

Differentiating ([^5)) with respect to the chemical po- 
tential and taking into account (|43p. we can verify the 
vahdity of (|44|. 

From (IMl) it is clear that the meaning of the bulk den- 
sity of the ri-potential can be assigned to P* . However, it 
will be seen from the following that the "local pressure" 
is not the only candidate for this role. 

Using (pS)) and taking into account ([TU| and (|M|) . it 
is easy to obtain the expression for the nonlinear surface 
coefficient 



+00 



[0P~P*{x)]dx, 



(96) 



which defines it in terms of 
similarity of dHU) and ((75)) . 



'local pressure". Note the 



3.10. "Mechanical definition" of 7 

By this definition is meant the procedure of compression 
and expansion of the system in two different directions 
with the total conservation of a volume and changes in 
definite internal surfaces [25l p. 43]. We apply this proce- 
dure in a modified form to the case of a wall/fluid inter- 
face. 
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Figure 1: The system for mechanical definition of the sur- 
face ri-potential of a fluid. The parameters of 
the lower surface differ from those of the remain- 
ing faces. 

We assume that the physical surface is inside the system 
at a microscopically large distance from its boundaries, 
and, at the same time, this distance is much smaller than 
the sizes of the system. 

The procedure is as follows (see FigH]). First, the upper 
plane of the system is shifted downward by 8h {6h > 0). 
The work done on the system at this stage is given by 



,51^1 = 



P{l-2ltf +AI 



/ Pstidx 



5h, 



(97) 



where Pst is the pressure in a narrow strip near the edges 
of the planqj, which differs from the pressure at the center 
of the plane (in the bulk). The subscript s indicates that 
this pressure, generally speaking, does not coincide with 
the tangential component of the pressure tensor. The fact 
is that this component is calculated for the effect exerted 
by the neighboring volumes of a fluid on each other, and 
in this case we deal with the impact of a fluid on the 
transverse wall. As this takes place, the whole picture at 
an atom level is highly distorted. 

At the second stage, the system is allowed to expand 
laterally 



and +00. Thus, it can be concluded that, for the fluid, 
the elementary work is expressed as 



SW = 5A 



f Pstdx + I {P- Pst)dx 



(103) 



5W2^ 



P{l-2k){h~2lt) 



It It 

+{l + 2h) I P,t,dx + l I Pst.dx 



(98) 



SI. 



Here, Pst2 is the tangential pressure acting along the 
bottom edge, and 61 > 0. 

The condition of conservation of the volume is not quite 
trivial: we must assume the constancy of the volume 
bounded by the planes of the arbitrary localization de- 
fined by the coordinate x' measured from the boundaries 
of the system 



where the tangential pressure Pst acts along the bounding 
planes. 

Since the process is conducted at constant N and T, the 
work should be equal to the change in the free energy F = 
F{V,T, N). Some complication is the pressure variation 
during the process, which is, generally speaking, due to 
adsorption. By varying 

F ^ -~PV + 7^ + /i(A^6 + Ns) (104) 

and taking into account that 5Ni, = —SNs, we obtain 

6F = -V5P + "fSA + ASj+{Nb + Ns)6fj,. (105) 

In the right-hand side of the equation all terms are can- 
celled in pairs, except the second term, and finally we have 



SF = -fSA 



(106) 



as it should be if 



(/ - 2x')^Sh ^{l- 2x'){h ~ 2x')dl, 
or, retaining the first-order terms, 

,1 — h\ 5h 



61 



l + 2x'- 



h 



(99) 



(100) 



-t-CXD 



Then, in the same approximation, we have 

6W = 6Wi + 6W2 



Ml 



■6A-i 



- f Psudx+ f{P-PstJdx 

an' 

x' It 

- I Pst,dx+ fiP-Pst,)dx 



(101) 



= - f Pstdx + f {P- Pst)dx. (107) 



Thus the "mechanical definition" of 7 also gives the 
same structure, which is already familiar to us (j78p . (IM)) . 

3.11. Tangential force and "local pressure" 

Let us show a consistency of the approaches ([M]) and 
(|107p . We differentiate (I104p with respect to the area 
using (P5|) for 7 



dA 



Ny 



-^ll+^ + ^ll + (^^ + ^^)|l- (^0«) 



where 



6A1 =\-2l + — \6h and M2 = —6h (102) 

\ h J h 

are changes in the areas of the first and the second types 
during the process, respectively. 

In this example, the shape of a parallelepiped was cho- 
sen not by chance: two parameters - / and h - enable 
Ai and A2 to be varied independently. This allows us to 
consider only one type of the area to be variable. 

By virtue of the assumptions made, the lower and the 
upper limits of integration can be considered equal to —00 



The first term in the right-hand side is cancelled with 
Niydfi/dA in the last term, and the third term is cancelled 
with Nsdfi/dA. This is easily obtained by going over to 
the derivative with respect to the chemical potential in 
the integrand and using (j43p . and we finally have 



dF_ 
'dA 



N.V 



P*(x)dx + [P- P*{x)]dx. (109) 



Thus, we can identify (|M|) and (I107p . and immediately 
obtain 



^For simplicity we assume that all the lengths of the transition re- 
gions lt,l3t,l* ,lu, connected with Pt,Pst,P*,S respectively, are 
identical, unless otherwise stated. Using the averaging integrals, 
it can be assumed that the integration is performed over the 
maximum of these lengths. 



Pstdx 



Pst = P* 



P*dx 



(110) 



(HI) 



where the bar denotes averaging over the transition re- 
gion. 
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In reality, in the absence of clear boundaries of the tran- 
sition layer, the accuracy of (jllip is determined by the 
proximity Put and P* to zero and P at the boundaries 
of the averaging interval. Nontrivial part of the equation 
is determined just by the difference between these quan- 
tities and their limits. Thus, the value of the averaging 
interval in (jllip is determined by a compromise between 
these two factors. 

3.12. Pressure tensor and near-wall distortions 

Let us now prove a similar equality for the tangential 
component of the pressure tensor 



Pt = P*. 



(112) 



We use the identity of the problem of inhomogeneity 
in the field with that of the liquid-gas interface for this 
component li\. Let us write the expression 

Pt{xi)^gikBT-^J{y2-yi)g^^^^dr2, (113) 

which is identical to the expression used in many papers 
[3j[T3l[T5l|2T] . Here Uij is the interaction energy between 
the i-th and the j-th particles, y is the coordinate parallel 
to the surface, and the field gradient, as before, is directed 
along X. 

A brief outline of the proof is as follows. 

Equations (|5^ and (J113I) are integrated over the coor- 
dinates of the first particle, and the integrands are then 
expanded in a series in the activity using (|20l) . Finally, 
their identity is proved by equating the coefficients at the 
same powers of z. 

A proof of the equality of the integrals of Pt and P* 
over volume is given in Appendix [C] 

We prove that in the region of homogeneity (no field) 
Pt — P* — P. Since all the directions in this case are 
equivalent, for a; and z we can write equations similar to 
(|113p . Composition of these three equations gives 



iPt = iQksT - 2 / £'12^ (^2 - ^i)"a^ 



, sdui2 , , .dui2 

(y2-yi)^ h(z2 -zi) 



dy2 



dz2 



dro 



(114) 



Collapsing the expression in brackets, we arrive at the 
well-known formula for the pressure 



Pt = gkBT 



\l 



(2) dui2 , n 

£'12^12^ dr2 =P 

dri2 



(115) 



H p.190]. 

For P* it is immediately obvious that for 9^ = \, i 
1,2,..., (112) reduces to ([TO]) or 



P* ^P 



(116) 



in the region of homogeneity. 

Thus, by virtue of equalities (|115p . (J116I) in the region 
of homogeneity and the equality of the integrals of Pt and 
P* over the volume of the system, we can conclude that 
pT^ is valid. 



As in the case of Pgt, the transition layer in (I112p is 
determined by the proximity to and 1 on the edges of 
the interval. 

There are some reasons to believe that the local equality 
is not valid, i.e., generally speaking, 



Pt 7^ P*, 



(117) 



but a rigorous proof of (|117p for these systems is beyond 
the scope of this paper. 

Also we shall not deal with the question of pressure 
tensor ambiguity, because only its zero moment is taken 
into consideration [2T1 . 

Thus, we have established the relationship between the 
tangential pressure and the behavior of near-wall distor- 
tions (oscillations) 

g{x)dx ^ jj- I Pt{x)dx (118) 



or 



g 



9/i ' 



(119) 



where the compliance with the macroscopic relation (J12p 
also takes place. 

Simultaneously we obtained the equality for the tan- 
gential force acting on the wall 



Pst — Pt, 

which is not obvious in advance. 

Finally, expression (j95p can be written as 



(120) 



X oo 

.M^-/p,Wdx + /[p-p.Ml.. ,121) 

— oo x' 

This expression extends the standard Kirkwood-BufF 
formula ((T|) to the case of the wall/fluid interface for an 
arbitrary position of the dividing surface. 

3.13. Equivalence of "adsorption" and "surface 
tension" approaches 



Equality pi2p allows us to prove the identity of the "ad- 
sorption" or "claster expansion" approach [fliSEIS on the 
one hand, and the "surface tension" or "pressure tensor" 
approach on the other hand |13ffT^[T7[|18[[^ . It suffices 
to show the equivalence of expressions ([55|) (with v in the 
form (EHl)) and (fnTj) . 

Canceling identical terms, we obtain the equality con- 
dition for these expressions in the form 



-|-CX) 

/ {0P - Pt)dx. 



(122) 



Comparing (I122p with (IM|) and taking into account 
(I112p . we can see that this equality is valid. 

In view of (j56p . we have another expression for the spe- 
cific surface f2-potential 



+ 00 

-pkh{x')+ f {ep-Pt)dx, 



(123) 
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which, of course, is identical to expression (I12ip but with 
the terms linear and nonlinear in pressure explicitly sep- 
arated from each other. Replacing Pt by P* in (J123I) . we 
return to the approach in terms of "adsorption" , replace- 
ment by Pst leads us to the "mechanical definition" of 
7- 

3.14. Near-surface virial expansion 



We can perform integration over dxi similarly to (|26p . 
Making the change of variables x'j^ = xi, r^ — ri — ri 
{i — 2, . . .n) and using the invariance oiU^ \. under trans- 
lations, we get 



dn = 



"2...n "0,2...n"'^2---«ni, 



(131) 



where 



Equation (|119p and properties (|112l) . (|120p allow us to 
construct near-surface analogs of the virial expansion. 

As a preliminary let's consider a case of a homogeneous 
environment. Then we have a system of equations (llOp . 
(jlip . which can be written as 



*2...n 



h 



/n 
O{x)\T[0{x + Xi) dx. (132) 

1=2 



(3P = 



n=2 

oo 

z -\- 2^ nbnZ^ 



where 



bn = — U["'}„dr2...drn. 



(124) 



(125) 



By virtue of the identity of the functional relation P 
and g in (I124p via z on the one hand, and P and 'g in 
(|129p via z0 on the other hand, we can immediately write 
down the analogue of (I126p 



I3P* 



k=l 



rSkt 



k+l 



(133) 



We emphasize that the integration in (I125P is performed 
over infinite space ^^. 
Eliminating z, we obtain 



where 
Sk = 



E(-i)^^--^^^^ 



^^y- TT ijdjY 



n 



rrii 



pP = g-J2- 



k 



-13kg' 



k+l 



(126) 



(134) 
and condition (I128p is preserved in its original form. 
Finally, taking into account (|112l) . (|120l) . we can write 



k=l 



where 



PPt^^- 



ijh 



■5-^ ^ 1 (fc ^ 1 + y^,- '71,)! T-r 

Pk^Y.i-^)^'"'' ^r^^n 

m 3 '' 

(127) 
are the so-called irreducible cluster integrals. The sum- 
mation is performed over all sets of a nonnegative rrij, 

satisfying 

fc+i 

Y^ij - l)m, = k (128) 

i=2 

H p.144]. 

As is easy to see the condition (|128p ensures the same 
dimensionality of different terms in the sum (|127p . 

In the presence of the field, we must start from the 
equations (1921) . (TTSl) . Averaging them on the transition 
layer, we obtain the system 



E' 

k=l 



k 



rSkQ' 



k+l 



and 



PP,t^g-J2- 



T^fc + l 



fc=l 



1 



SkT 



(135) 



(136) 



Expression (|135p and (I136P is preferable to (jl33p : a Pt 
and Pst can be extracted from computer experiments. In 
addition, Pt can be calculated, because it has a represen- 
tation in quadratures (J113p . 

From the point of view of a nontriviality, averaging in- 
terval in (I135P , (|136p must be chosen minimal, but so that 
the equalities (lllip . (I112p are still valid. 

Thus, the connection Sk with dj defined by p34p . For 
example, the first few relations are 



/3P* = ze + J2dnize)" 

n=2 

oo 

'g = Z0 + Yndnjzey 



Si - 2d2 




S2 = 3^3 - 


-6d^ 


(53 = 4^4 - 


- 24(^2^3 



(137) 



80 



(129) 



where 



dn = „ / / I J|6'JzY{"|'„dxidr2...dr„ 

1 1 • if \7 1/ 1/ ■ -1 



(130) 



Here the integration on the coordinate of the 1-th par- 
ticle is performed over the transition layer, and on the 
coordinates of the other particles - over the infinite space. 



and they are identical to the corresponding relations be- 
tween Pk and bj [H p. 144]. 

4. Discussion 

4.1. Localization of the dividing surface 

As is seen from the foregoing, many difficulties arising 
in studies of adsorption are related to the fact that the 
authors combined two surfaces: system and field. Their 
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separation made in the present contribution by transfer- 
ring the field into the system immediately clarifies the 
problem. 

It is variant (l3T|) which is used in papers |4l p. 424], 
and [21 section 33] , as well as in earlier works given therein. 
This follows from the fact that the derivation procedure 
employed in these papers assigns the surface part to the 
liquid volume. However, this was done incorrectly. 

It turns out that calculation of surface terms calls for 
integrating beyond the limits of the system. The reason 
is that surface effects are determined by the difference 
between inhomogeneous and homogeneous number densi- 
ties. As for homogeneous number density, it requires that 
integration be made over correlation volume to both sides 
symmetric about the surface. 

In the above works integration was rigorously restricted 
to the system volume, i.e., it did not go beyond the limits 
of its surface, and therefore could not give a correct result. 
(Note that this is valid for all terms of the series except 
the first term). 

This approach is also used in current works [T? . 

This question was discussed in the literature, and the 
above incorrectness was eliminated [5]. However, the form 
of the first term remained the same. The dividing surface 
was fixed at the first atomic layer of a solid: the employed 
method is unsuitable for other variants [5]. 

Even in the case of correct calculation this variant is 
inconvenient, since it assigns surface particles to the liq- 
uid volume. In particular, given the attraction, i.e., real 
adsorption, in the first order in density the number of 
surface particles 



Ns^ gA I [exp(-/3w) - l]dx 



(138) 



may become negative for definite forms of the interaction 
potential, which is seen from (|138p and (03]). 

It can easily be understood, indeed, actually in expres- 
sions (PT|) , (I138P we assign the entire transition region to a 
homogeneous liquid, thus correction should be subtracted 
from the number of surface particles. 

Variant ([50)1 is preferable in the presence of real adsorp- 
tion, since it agrees with the common notion of particles 
adsorbed on a solid surface. In the first order in g it always 
gives a positive value of the number of surface particles. 
As follows from (l40ll. in this case 



(139) 



i.e., always positive. 

Here particles are divided into two parts: particles be- 
longing to the homogeneous liquid and those belonging 
to the transition layer; they are not taken into account 
twice as in (jl38|) . In other words, we add them in (I139p 
as external relative to the liquid, and in (I138P they are 
subtracted as internal. 

The attempt to directly identify these variants leads to 
two definitions of adsorption, "Gibbs" and "mechanic" , as 
well as to the introduction of the "effective potential" [5] . 



Ns ^ gA i ex.p{—/3v)dx, 



Variant (j32[) is inapplicable in the presence of real 
adsorption, since it assigns molecules adsorbed on the 
surface to the liquid volume, i.e., it is obviously a 
workaround, artifice in this case. Indeed, the presence 
of pronounced negative potential near the surface, as is 
readily seen from (l54l) , shifts the surface position xq deep 
into the solid. 

This variant is good when the potential of the surface is 
close to that of a hard solid or, speaking more generally, 
in the absence of a pronounced adsorption. 

4.2. Surface number density 

At arbitrary number density when the dividing surface is 
placed on external boundary of the field, expression (j43p 
takes the form 

CO 

gs = g{x)dx + / [g{x) — g] dx, (140) 

where the external boundary of the field in the second 
integral is placed at the origin of coordinates. 

Here the two terms are defined most sharply. The first 
term in the right-hand side of (|140p corresponds to the ad- 
sorption in a classical understanding as particles residing 
on the surface of a solid, is linear in the number density, 
and is defined by the region of direct action of the wall 
field on fluid particles. 

The second term is quadratic in the number density 
(which follows from P^ ). and is specifled by indirect ac- 
tion of the field by its transfer through the layer of ad- 
sorbed particles. It is related to interparticle interaction 
and geometric displacement of particles from the region 
of "dead" volume near the surface. 

As is should be, for the hard solid only the second term 
works in (jl40p . and it reduces to (|^5)) . while for macro- 
scopically smooth field - only the first term, and (jl40p 
transforms into 



g{x)dx, 



(141) 



which, as is easily seen, agrees with ([77)) . 

The second term in (|140p may be called "indirect ad- 
sorption" unlike ordinary one (the first term) that may 
be called "direct adsorption". Though, these two terms 
are not strictly separated. 

Some remarks concerning the dividing surface. Differ- 
entiating expression ([65)) with respect to the chemical po- 
tential, one can see that in this case the most natural 
choice is the Henry constant of adsorption is equal to zero. 
Exactly this choice of the dividing surface is used in [S] . 
Recall that it corresponds to the inclusion of near-surface 
volume AR, inaccessible to fluid particles, in the solid 
volume. 

On the contrary, in a great number of computational 
and other works it is common practice to assign "dead" 
volume to a fluid one, which results, as is evident from 
([7^ . in additional negative term —R in the above- 
mentioned constant yjj . Note that results of paper [11] 
for surface number density (both theoretical and those 
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obtained by computer simulation) agree with (|74l) . This 
is seen from the comparison between the derivative of this 
expression with respect to chemical potential and the ex- 
pansion of two variants of PercusYevick equation into a 
series in density [TT:. 

In conclusion note that unlike the case of the liquid-gas 
equilibrium system, in the variant at hand imposing the 
condition Ng — 0, as is seen from (l43l) . results in that the 
position of the dividing surface becomes dependent on the 
number density. Probably, in the given case there is no 
physical meaning of such a separation. Instead, in some 
situations one can use condition (j33p . 

4.3. Tangential force vs pressure tensor 

As was mentioned above, we cannot identify the specific 
tangential force and tangential component of the pressure 
tensor just from general considerations. In other words, 
we cannot a priori believe that Pgt = Pt and even that 
Pst — Pt ■ There are two reasons for it. 

First, Pt is introduced for the conditions, which ensure 
a uniformity along the surface. It specifies the interaction 
of the neighboring volumes of the fluid. On the contrary, 
entering of Pst implies the presence of a transverse wall, 
i.e., the involvement of a corner structure. It is clear, that 
in this case particles distribution has nothing in common 
with that in the previous case. 

Second, Pt is not the force, but the time derivative of 
momentum [Tl] . It consist of the force and the kinematic 
terms. The special theorem (or theorems) is needed to 
certify that this momentum flow is equal to the corre- 
sponding part of the force. 

The proof that P st — Pt given in the present paper is 
one of such theorems. 

4.4. Adsorption and surface tension 



Obviously expression (|44)) is of a universal nature. So the 
surface number density (which just specifies an adsorption 
in a broad sense) is closely connected with the specific 
surface fi-potential 7. 

As follows from the "mechanical definition" of 7 (Sec- 
tion I3.10p and theorem (I120p , the surface forces are de- 
termined by the same value. 

Thus, we may state that the surface tension (also in 
a broad sense) for the case in point is closely connected 
with adsorption. 

Therefore, the proof of the identity of the two ap- 
proaches denoted as "adsorption" and "surface tension" 
performed in the present contribution is primarily related 
not to these notions but to methodically different direc- 
tions of the research. The first approach initially appeals 
to the systems in the equilibrium, while the second one 
starts from the analysis of the force dynamics |13j . 

On the other hand, the reasonability of using the terms 
"adsorption" and "surface tension" in the given context 
is undoubted, due to the features of the evolution of the 
above-mentioned approaches. 

Note that "direct" and "indirect" adsorptions are not 
strictly distinguished, since both terms are involved in 



the tangential pressure. Besides, as is seen from the com- 
parison between (P5)) and (|57|) . both of them, generally 
speaking, give a contribution to the nonlinear surface co- 
efficient V. 

4.5. Surface thermodynamic potentials 

Equation ([M| coincides with the expression for $7- 
potential of the liquid/gas system [3] ■ It is this potential 
which is most natural for the phase boundary. However, 
for the liquid/gas interface we can always choose the posi- 
tion of the dividing surface such that the surface number 
density goes to zero. In this case the surface f] - potential 
coincides with the corresponding free energy |3|. 

As already mentioned, in the case under consideration 
the surface number density does not go to zero, so 7(2:') 
is no longer a specific surface free energy. Indeed, 



F = n+fiN = -P{V-V')+-f{x')A+fi{Nb+Ns) = Fb+Fs, 
where 



Fb = -P{V - V) + ^iNb 



(142) 
(143) 



and 



Fs = A[-f{x') + ^iQs] (144) 

are volume and surface free energies, respectively. 

5. Summary 

1 . The general expression of the fi-potential (grand po- 
tential) of the system with locally introduced force 
field is derived (j25p . This technique makes it possi- 
ble to vary the position of the dividing surface (P^ . 

2. Volume and surface terms of the J7-potential cannot 
be unambiguously separated, and their interpretation 
is specified by the form of the solid-fiuid interaction 
potential pi)) - (15^ . 

3. In particular, depending on the situation, the surface 
term can begin both with linear and with quadratic 
term in the activity ([50)) . (|32p . which corresponds 
to the presence or to the absence of adsorption in 
ordinary understanding. 

4. Solid phase contribution to the surface number den- 
sity is defined by (?(x), and fluid contribution - by the 
difference q{x) — g (|^ . 

5. The general expression for the specific surface J7- 
potential 7, composed of the terms linear and non- 
linear in pressure is obtained (j56p . The linear term is 
determined by the product of pressure and the Henry 
adsorption constant (|50)l . and depends on the loca- 
tion of the dividing surface. The nonlinear term (the 
nonlinear surface coefficient - v) has several universal 
representations dH]), (H?]), dMl, ([T^ - 

6. The surface number density, similarly to the 7, 
breaks up into the term linear in the bulk density 
depending on the dividing surface position, and uni- 
versal nonlinear term (157 



14 



7. The Henry constant of adsorption depends on the Appendix A Factors 

position of the dividing surface, and is determined 
by the universal expression x' — xq (|53p . 



A.l Ursell factors 



?. The expression for the nonhnear surface coefficient v 
is obtained, in the form of a series in powers of the 
activity - the analog of the Mayer series for pressure 
((27)) . The coefficients of the series are given as in- 
tegrals of the products of Ursell factors and simple 
multipliers, that depend on external potential 



9. For the hard solid case (section 13.51) the coefficients 
of the series are the first moments of Ursell factors 
taken over certain regions ([7^ . 

10. For macroscopically smooth field of arbitrary config- 
uration (section I3.6p the expression obtained for v 
(j76p agrees with thermodynamically reasonable ex- 
pression for rj-potential ([77]) . 

11. Probably, the system with exponential wall- fluid po- 
tential, can be used as a convenient model for the 
calculation of surface terms (section 13. 7p . 

12. The expressions for v and fi-potential of microscop- 
ically inhomogeneous held of periodic surface struc- 
ture are deduced fsection l?^ . 

13. The "surface tension" approach produces the expres- 
sion for V in quadratures (jl22p . which permits to op- 
erate with 7 as a whole (|123p . 

14. The "adsorption" approach produces the expression 
for V in quadratures, when "local pressure" (j96p is 
used. 

15. On average over the transition region near the surface 
the tangential component of pressure tensor (J113p co- 
incides with the pressure exerted on the transverse 
wall (HZl) and with the "local pressure" ([92]). This 
fact together with arbitrary position of the separat- 
ing surface allows one to consider new versions of the 
Kirkwood-BuflF formula dSS]), (fTUT)) . dHH). 

16. The ri-potential of a statistical system in a force field 
of arbitrary configuration is defined by the integral of 
"local pressure" over the entire volume ((89)) . (l92|) . In 
view of the previous item, it is also defined by anal- 
ogous integrals of Pt and Pst, when these quantities 
are valid. 

17. On average over the transition region the relation 
between the near-surface number density and tan- 
gential component of pressure tensor corresponds to 
ordinary macroscopic relation between number den- 
sity and pressure (|119|) . 

18. For the problems under study, "adsorption" and 
"surface tension" that evolved as independent lines 
of investigations are completely equivalent in their 
common domain of existence (section 13. 13p . 



19. The near-surface virial expansion (section I3.14p de- 
termines the equation of state of near boundary "two- 
dimensional" fluid. 



Ursell factors may be deflned by the equality: 



{n} a=\ 

I 

1 < ka,<k, J^ka^k, exp(-/3C/i) = 1, (A.l) 

a=l 

where by {n} we denote some partition of the given set 
of k particles with the coordinates ri, ...r^ into disjoint 
groups Tia, I is the quantity of groups of a particular par- 
tition, ka is the size of the group with the number a. The 
sum is taken over all possible partitions, and the meaning 
of the condition exp(— /3[/^) = 1 is evident: unit groups 
in the given case make no contribution to the products. 



For example, the first several U[ \. are 



U^^^ 


= 1 


,/(2) 
"l,2 


= exp(-/?C/i%) - 1 




= exp(-/3C/32 3) - exp(-/3t/2 




- eM-PUli) - exp(-/?C/2'3 



(A.2) 



As is well-known, one of the generating functions for 
Ursell factors is the logarithm [55]. However, the division 
of the series in ^ for the case A: = 1 by Sy in form ([7]) 
shows that there exists the second generating function 
- fractional - in the form of the series ratio ([SJ. This 
follows from the familiar expansion of the number density 
for homogeneous system pip that is valid far from the 
system boundaries. 

In Appendix IB. II the existence of fractional generating 
function will be proved. 

A.2 Partial localization factors 

As far as we know, these functions were introduced for 
the first time in [?7] . 

Some particles involved in these functions do not cause 
the decay when moving away (the delocalized group) , and 
some of them - do (the localized one). 

Introduce the notation 



B\ 



{in^k) 



-k^ 



(A.3) 



where the superscripts m and k define the quantity of 
delocalized and localized particles, respectively {m — 
1, 2, 3, . . . , fc = 0, 1, 2, . . . ). The subscripts denote the co- 
ordinates of particles, with the first m particles being con- 
sidered delocalized, and the rest - localized. 

These functions are similar in structure to the Ursell 
factors of fc -|- 1-th rank, with the proviso that in the con- 
struction by the type of (jA.ip the first m particles (de- 
localized) are treated as a single compound particle. In 
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(m,k) 



other words, define B\ 'J^j^^^ by the equahty 



n(m,k) 
'-^l...m+k 



{n} 



(A.4) 



-rive at 






"l...m+fc 


_ „(m+fc,0) 
- "l...m+fe 


fc 
\ ^ \ ^ »3(n,0)t,(m,fe-n) 

/ . / . "l...n»n+l...m+fc' 
n— 1 samp 



n exp[-/3C/ 



fe„+(im-l)5aT 



{no 



a=l 



l<fca<A: + l; ^/ca = fc + l; exp(-/3t/^) = 1, 



(B.l) 



where m > 1, fc > 1, and the internal sum is taken over 
samplings of localized particles only {n from fc). Relation 
(jB.ip is proved either by direct exhaustion of partitions in 
accordance with (|A.4p . or by repeatedly substituting the 
expression for B[^'^^i, in the right-hand side of (jB.ip . 



where designations are analogous to (jA.ip on condition 
that summation is taken over all possible partitions of 
the set of fc + 1 particles among which one particle is com- 
pound. Sar is the Kronecker delta, and r is the number 
of the group involving a compound particle. 

The generating function for ^Si „ i_j. is the distribution 

function QqI^ ^ of GCE type. Expanding the partition 
function Sy in ^ and dividing the series, we obtain (ITl)) . 
The proof of this relation is given in Appendix IB. II 



Paper [55] gives another recurrence relation for Ursell 



functions which in terms of B 



{rn,k) 



.m+k 



looks like 



"l...fc+l 



„(fc+l,0) 

"l...fe+l 



k 

E 

n=l 



„(»i,0)„(l,fe 
"l...n"Ti+l 



+ l...k+l 



^I'm+k th^t ^""^ the first in localized group: 



B]; 



(m,0) 



exp(-/3t/i" 



(A.5) 



(m4) 



- exp(-/3t/r.+i,„ , 2) - exp(-/3t/i".„J 



X exp(— /3L/, 



m+l,m+2) 



2exp(-^C/™ 



perm 

(B.2) 
where fc > 1, and square brackets denote averaging over 
permutations of particles. 

In principle, this is the same equation (IB.ip at ttt, = 1 
but written in a symmetric form. Note that in the given 
case there is no necessity in symmetrization - equation 
(|B.ip holds rigorously in asymmetric form as well. Never- 
theless, (|B.1I) can also be brought into a symmetric form, 
however, this calls for averaging over permutations of both 
the sums and the left-hand side of the equation. 

B.2 BBGKI equation 

In the presence of external field this set of linking equa- 
tions may be written as 01 p. 205] 



and in delocalized one: 



„(i,fc-i) 

"l...k 



U 



(fc) 
l...fc' 



^lSl...rr. 



(A.6) 






iMl, 



including, for homogeneous medium 



-/5 



Qi^^ll"^ lUlm+ldr^n+l- (B.3) 



B\ 



(1,0) 



1. 



(A 7) This relation, after substituting expressions (|20p in it 
and equating coefficients at equal powers of z, generates 
It is obvious from (|X5|) . and (|X6)) that partial localiza- the differential recurrence relation 
tion factors generalize the concepts of Boltzmann factors m 

and Ursell factors, including them as the limiting cases. ViB^'^,^ ~ — /3;BJ™^" , „y^ ViMii 

For factors Bl"''^,-^. a number of recurrence relations 
hold, which ensure the existence of various physical links; 
some of them are given in Appendix iBl h /_^ ^ii..m+n 

2— m-t-l 



ra^n 



i=2 
j(m+l,ri— 1), 



-/3 E M™rCVi.n, (B.4) 



Appendix B 



Recurrence relations for 

(m,fc) 
m+k 



B\ 



Many relations and operations of statistical mechanics 
are provided by definite classes of recurrence relations for 
^im+k- ^^^ brevity, we shall say that an operation gen- 
erates a recurrence relation or a class. In this contribution 
only some of them will be considered. More detailed in- 
formation can be found in |22) . 

B.l Correspondence to the definition 

Equating expressions dU and ([T4)). expanding the series 
for Sy, and performing multiplication of the series, we 



which we use to prove (|112l) . (jllSp . In the second term 
on the right all permutations of particles of the localized 
group (n) with the second particle are exhausted. 

Appendix C Proof of (Pi) = (F*) 

To prove this equality, consider the system inside the 
volume shaped as rectangular parallelepiped strongly 
stretched along the y axis (see FigH]). The system includes 
the interface plane, and one of the system boundaries lies 
fairly deep inside the field region. Wall/fluid interaction 
is still defined by the functions 6. 
Let us integrate expressions 



((92|) and (|TT3)) with respect 
to the coordinates of the first particle. The integration 
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g(i,n) g^j,g exactly the factors of this kind, the integral 
term goes to zero. Thus, 



Figure 2: The system used to prove {Pt) = (P*)- Fluid is 
at the top of the interface. 



is performed over the whole system volume denoted by 
dashed line in Fig. [21 We put that this volume is defined 
by the characteristic function -01. 

Expand the integrands into a series in the activity using 
([20)). Equating the coefficients at equal powers of z, we 
have 



keT f 01 




^l'.''^k+2^'^i--d-'^k+2 


(fc + 2) 
2 


/^l 


fc + 2 n 

i— 1 J 


, NK(2,fe) dui2 



(C.l) 



dri...drk+2, 



where k = 0, 1, ... . Here we equated the coefficients at 
z'^"'"^. So, it is necessary to prove this equality. 

The kind of the recurrence relation ()B.4I) required to 
prove (|C.ip is of the form 



^Bl gy/i...fc+2 



fc+2 

/ . "21. .fc+2 
i=l,i^2 



?(2,fc) 



du2i 

dy2 



(C.2) 



Multiplying ([C.2I) by {y2- yi)^idi---0k+2 and integrat- 
ing over all coordinates, we obtain using the symmetry un- 
der permutations inside localized and delocalized groups 

of S(™:») 



ksT / -01 



fc+2 



iy2 - yi)—B[^'^^^^^]dri...drk+2 



= Ii+ kh + kh 



(C.3) 



where 
h = 

h = 



i'l 



^1 



fe+2 



(2/2 - yi)B\^^'f,^^^^—dri...drk+2, 



fe+2 



j(2,fc) 



9^2 



C?M21 



(C.4) 



(2/2 - 2/3)^2,1... fc+2 -^ — dri...drk+2, 

(C.5) 



-^3 == / (03-01) 



fc+2 



/ \f2(2.fc) ^"21 , , 

(2/2-2/3)^2,1. fc+2 ^^^C^'^l-f^'^. 



.'=+2 ay. 



fe+2- 



(C.6) 
Here in the integrals I2 and I^ the change of variables 
I'l ^ ^3 is performed. 



The integral in the left-hand side of equality (jC.3|) may 
be taken by parts, since here factors 0i depend solely on 
Xi. Wherein, due to locality of Ursell factors, and factors 



- ksT / 01 



fe+2 n 



(2/2 - 2/i)^si'.l+2 dr-i...drfc+2 



= ksT / 01 
It is readily seen that 



fe+2 



j(i,fe+i) 



^l'.fe+2*l-'^^fc+2- (C.7) 



/i=2/2. (C.8) 

For this purpose, calculate the difference 

Ii-l2= ipi Y[0^ {y3-yi)B^^f'^j^^2-p^dri...drk+2- 

(C.9) 
Making the change of variables ri ++ r2, we see that 



the right-hand side of <\C.9\i coincides with I2, so we have 



Finally, consider the integral I3. Obviously, the factor 
(03 — 01 ) is nonzero only in the case where the first and 
the third particles are on different sides of the system 
boundary. It is seen better from the expression 



03 - 01 = 03X1 - 01X3, 



(CIO) 



where we used the identity 0^ -I- x^ = 1. 

When the 1-st or the 3-rd particle moves away from 
the boundary of the system, the integrand decays rapidly 

(2 k) 

owing to local character of the product S, 1 fc_|_25M2i/52/2. 
Hence, the integral is proportional to the surface of the 
parallelepiped at hand, but not all of its faces make a 
contribution to I3. 

Consider the top face of the system, where due to re- 
moteness from the interface we can consider that 0i — 
l,i = 1,2,3.... Making the substitution (JC.IOI) . we ar- 
rive at the conclusion that condition I^ — demands in 
this case that 



?(2,fc) 



du2i 



03Xi(2/2 - 2/3)^2.1. ..fc+2^^; — dri...drk+2 



dy2 



(C.ll) 



?(2,fc) 



du2i 



01X3(2/2 - 2/3)^2,1. ..fc+2-^ — dri...drk+2 



However, obviously, equation ([Clip is true by symme- 
try when the y axis is parallel to the plane defined by 
the function 0. (To prove this it suffices to consider the 
mirror refiection operation with its mirror-plane coincid- 
ing with the upper face of the system.) Thus, we may 
conclude that contribution to I3 from the whole region 
of homogeneity above the interface equal zero except the 
side faces y = const. In the region of homogeneity be- 
low the interface this contribution goes to zero owing to 
9, -^ 0. 

So, only the side faces y = const and two narrow strips 
on the faces z = const make contributions to 73. The 
contribution of these strips is of higher order of smallness 
and may be neglected. Thus, we have the estimation 



Jzv 



(C.12) 
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The contribution of I3 tends to zero with the system tion in the form of (|113p is not accidental. We could omit 

length (along the y axis) tending to infinity. Here Szy is the term yi in the integrand, but it will lead to poor lo- 

the area of the face zy. calization of functions and complicate the proof. 

In view of (jC.Sp and (|C.12p . equation (|C.3p gives (jC.ip . In a similar way, but easier, it may be proved that rela- 

Note, that the selection of pressure tensor representa- tion dB corresponds to expansions (HU]) and (HIl), dH])- 
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